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Abstract 

The next-to-leading order inclusive cross section for large-px photo-production of charm 
quarks at HERA is calculated in two different approaches. In the first approach the charm 
quarks are treated as massive objects which are strictly external to the proton and the photon 
while in the second approach the charm mass is neglected and the c quark is assumed to be one 
of the active flavours in the proton and photon structure functions. We present single-inclusive 
distributions in transverse momentum and rapidity including direct and resolved photons. The 
cross section in the massless approach is found to be significantly larger than in the massive 
scheme. The deviation originates from several contributions which are disentangled. We argue 
that large-p^ photoproduction of charm quarks at HERA will be sensitive to the charm content 
of the photon structure function. 
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1. Introduction 



Charm production in high energy ep collisions at HERA is dominated by photoproduction events 
where the electron is scattered by a small angle producing photons of almost zero vir tuality (Q 2 ~ 0). 
In leading order (LO) QCD the main process is photon-gluon fusion where the photon interacts 
directly with the gluon from the proton producing a cc pair in the final state (7 + 5 — > c+c). Besides 
the direct photoproduction channel charm production at HERA can proceed also via the resolved 
photoproduction process where the photon behaves as a source of partons which interact with the 
partons in the proton, as for example in g + g — ► c + c. It is known that for the photoproduction 
of light quarks q (i.e. q = u,d and s) and gluons the cross section of the resolved process is much 
larger than the cross section for direct photoproduction Only at rather large transverse 

momenta, p±_ J> 30 GeV, the cross sections for the two types of processes are comparable. This 
has been confirmed in various ways by measurements of the ZEUS || and HI || collaborations 
at HERA. Actually a recent comparison of the ZEUS experimental data for the inclusive one- 
jet cross section JjJ with the next-to-leading order (NLO) predictions superimposing resolved and 
direct photoproduction shows satisfactory agreement in absolute normalization, p± and rapidity 
dependences Q. These predictions are meant to be valid only for large enough p±, i.e. pj_ J> 5 GeV, 
so that soft processes, which can not be calculated in perturbative QCD, are negligible. In the 
calculations the authors assume four active flavours, q = u,d,s and c, which are taken to be massless. 
This means that also charm quarks are produced and counted as jets. Of course, the c quark is 
also an ingoing parton originating either from the proton or the photon in the case of the resolved 
contribution. To neglect the charm mass m in these calculations is a reasonable approximation since 
in large-pj_ jet production the effective scale is p± which is much larger than m — 1.5 GeV. The 
same approach is followed for predicting inclusive cross sections for the photoproduction of light 
hadrons, like tt's, K's etc., when the produced massless c quarks fragment into the light hadrons 
similarly as the produced u, d, s quarks and the gluons Q. In all of these theoretical investigations 
the c quark is one of the massless active flavours in the photon and proton structure functions. In 
the following we shall refer to this as the massless charm scheme, which we assume to be valid in the 
region of large transverse momenta p± 3> m. In this scheme the small pj_ region is not calculable. 
The cross section diverges in the limit pj_ — > and total production rates can not be predicted. 



In the massive charm scheme which has been adopted by many authors (for a review see |T( 
the charm mass, m 3> Aqcdj acts as a cutoff and sets the scale for the perturbative calculations. 
Similarly to the case of large-pj_ jet production the cross section factorizes into a partonic hard scat- 
tering cross section multiplied by light quark and gluon densities [11]. Inherent in this factorization 
is the notion that the only quarks in the hadron and the photon are the light ones. Thus in the 
massive charm scheme the number of active flavours in the initial state is equal to rif = 3 and the 
massive c quark appears only in the final state. This approach has the advantage that not only 
various distributions, like in rapidity and/or transverse momentum, can be predicted but also the 
total cross section. In LO direct production is described by the partonic reaction 7 + g — > c + c 
while the resolved contribution involves the channels q + q — > c + c and g + g — > c + c, where q(q) are 
light (massless) quarks. NLO corrections have been calculated to these processes and found to be 
substantial p^Jl3[| . A comparison between the theoretical NLO results and experimental data on 
photoproduction of charm quarks in low-energy jp and 77 collisions has shown reasonable agree- 
ment [14]. Recently, the total charm production cross section at HERA has been measured [15] and 
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found to be well described by the prediction of Frixione et al. |l6|]. Depending on the choice of the 
photon structure function the resolved cross section seems to contribute only a fraction to the total 
charm production cross section |l6[ . 

One might expect that the massive approach is reasonable only in those kinematical regions 
where the mass m and any other energy scale like p± are approximately of the same magnitude and 
significantly larger than Aqcd- Under these circumstances the charm mass is used as the scale of 
a s and the quark and gluon densities of the photon and the proton. In next-to-leading order terms 
~ a s (fi 2 ) ln(p^_/m 2 ) arise from collinear emission of gluons by a heavy quark at large transverse 
momentum or from almost collinear branching of gluons or photons into heavy quark pairs. These 
terms are not expected to affect the total production rates, but they might spoil the convergence 
of the perturbation series and cause large scale dependences of the NLO result at p± 3> m. In the 
massive approach the prediction of differential cross sections is thus limited to a rather small range of 
p± ~ m. The proper procedure for p± ^> m is to absorb the terms proportional to ln(p^/m 2 ) into 
the charm distribution functions of the incoming photon and proton and into the fragmentation 
function of c quarks into charmed hadrons. Of course, to do this absorption one needs a charm 
contribution in the structure function in the first place. An alternative way of making predictions 
at large p± is to treat the charm quarks as massless partons from the start. The mass singularities 
of the form ln^^/m 2 ) are then absorbed into structure and fragmentation functions in the same 
way as for the light u, d, s quarks. We expect this massless approach to be better suited for the 
calculation of the differential p± distributions up to NLO in the region p± ^> m. Then the problem 
arises how to proceed in the intermediate range where p± > m. 

In order to investigate the region where p± > m we have calculated the differential cross section 
d 2 cr/dydp 2 L as a function of p± with the rapidity y integrated over the region \y\ < 1.5. We shall 
compare the results in the two approaches: (a) the massive charm approach with m = 1.5 GeV, 
in which we have computed the cross section for open charm production and (b) the massless ap- 
proach, where we have evaluated the same differential cross sections for inclusive charmed particle 
production. In both calculations we include direct and resolved processes and go up to NLO where 
necessary. The massive calculation is based on the work presented in [13] while the massless predic- 
tions use the results obtained in []l7j]. In the second approach we need the fragmentation function 
of the c quark into charmed hadrons. This is approximated by 5(1 — z) where z = pd/Pc is the 
scaled momentum of the charmed hadron, meson or baryon. With this choice the LO results in the 
massive scheme approach the LO massless results in the limit m — > if we restrict ourselves to the 
same parton subprocesses. They must differ, however, in NLO where the limit m — > is not possible 
due to the unabsorbed mass singular ln(p^_/m 2 ) terms. Of course it is no problem to incorporate 
more realistic c fragmentation functions in both schemes but the choice above should be sufficient 
to see the essential differences of the two approaches. A similar study of the production of large-p^ 
hadrons containing bottom quarks in pp collisions has been performed by Cacciari and Greco IS]. 
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2. Comparison of Results 



We have calculated the cross section d 2 a/dy dp\ in the HERA laboratory frame with E p = 
820 GeV protons colliding with E e = 26.7 GeV electrons travelling in the (+z) direction {y > 0). The 
virtual photon spectrum is described in the Weizsacker- Williams approximation with the formula 



f 7 /e(x) = £ 



i + d-^^ + a mlx ( i 




(i) 
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with Qmax = 4 GeV 2 , (5 min = m 2 x 2 /(I — x) and x = E^/E e in the interval 0.15 < x < 0.86 as in 



[15] which corresponds to jp cm. energies of 115 GeV < W < 275 GeV, m e is the electron mass. 
Since the average Q 2 ~ 10 -4 GeV 2 is very small, the photons are essentially on-shell, so that the Q 2 
dependence of a lp and the longitudinal contribution in (|l|) can be neglected. For the proton structure 
function we use the MRS(G) set IS] which describes well the proton structure at small x and is 
adjusted in the intermediate x range. In Fig. la the results for the cross section d 2 a/dydp 2 L averaged 
over the rapidity range \y\ < 1.5 are presented for the direct contribution in the massless and in 
the massive scheme. Both LO and NLO predictions are shown. In the LO cross section the same 
NLO parton structure function MRS(G) has been adopted, only the parton-parton scattering cross 
sections are evaluated in LO. The two-loop formula for a s is used with the A value taken from the 
MRS(G) fit: A^ = 255 MeV. Both NLO calculations have been performed in the MS factorization 



and renormalization schemes. The corresponding scales have been set to fi = \Jp\_ + m 2 . It is clear 
that in the massive scheme only three flavours are active in the initial state and in the evaluation 
of a s whereas in the massless scheme also the charm distribution in the proton contributes to c(c) 
production and a s is calculated using four active flavours. We observe that the predictions in the 
massless and massive scheme are very similar. In LO there is some difference, the LO massless 
cross section is approximately 20% larger than the massive cross section over the whole range of 
p± between 3 and 15 GeV. This small difference originates from the additional charm contribution 
in the proton and the value of a s . In NLO the direct contributions of the massless and massive 
theories yield almost identical results for the p± distribution in the range 3 < p± < 15 GeV. We 
have checked that this is true also for higher values of the transverse momentum up to p± ~ 30 GeV. 

The situation is completely changed for the resolved part of the c/c production. This contribution 
has been calculated with the NLO photon structure function of GRV |2(| in the MS scheme. The 
results for the p± distribution are plotted in Fig. lb. We observe that the NLO corrections in the 
massless scheme are large and increase the cross section by roughly 100%. In the massive scheme 
we have calculated only the LO result, which is between one and two orders of magnitude smaller 
than the LO massless cross section in the p± range between 3 and 15 GevS NLO QCD corrections 



do not change this order-of-magnitude suppression [16] and will therefore not be taken into account 
for the resolved contribution in the massive scheme. We have examined that the prediction for the 
massless resolved cross section does not depend strongly on the particular choice of the NLO GRV 
parametrization of the c distribution in the photon. The NLO GRV parametrization in the DIS 
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schemed yields similar results as well as the parametrization by Aurenche et al. [22]. 



1 Similar results have been obtained in the LO analysis presented by R. Godbole at the Photon 95 conference, 
Sheffield, UK, April 1995 @. 

2 To be consistent the appropriate changes have been made in the NLO direct part from the MS to the DIS 7 scheme 
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We may ask where the strong contribution of the resolved cross section in the massless scheme 
is hidden in the massive approach. It must be contained in the ln(p^_/m 2 ) terms in the NLO direct 
cross section. In the NLO massless direct cross section all these terms are absorbed either in the 
structure function of the proton, the photon or the fragmentation function. Therefore at NLO the 
definition of direct/resolved photoproduction is ambiguous and only the sum of both contributions 
should be considered. In Fig.lc we have thus plotted the sum of the NLO massive direct and the LO 
massive resolved contribution from Figs. la and b (denoted NLO/LO massive) and the same sum 
for the massless theory taken also from Figs. la and b (called NLO/LO massless). Considering the 
same order in a s these are the two predictions in the massive and massless theories which should 
be compared. We observe that the NLO/LO massless result is approximately 60% larger than the 
NLO/LO massive prediction. We have investigated this difference in some detail. To begin with 
one has to remember that in the massless theory the collinear singular contributions in the final 
charm states are absorbed into the fragmentation function. This subtraction is not performed in 
the massive theory. If we subtract the corresponding \n(p\/m 2 ) terms in the massive theory the 
NLO/LO massive prediction is increased by 25%, so that we are left still with a difference of 35%. 
This difference originates mainly from the charm content of the photon which does not correspond 
exactly to the ln(j?^/m 2 ) terms in the massive theory arising from the almost collinear branching 
of photons into charm quark pairs. We have examined this further by substituting the GRV charm 
photon structure function by the massive LO "Bethe-Heitler" approximation This change 

reduces the difference between the massless and the massive NLO/LO predictions to a level of 10%. 
Part of the difference between the two approaches can thus be attributed to theoretical uncertainties 
concerning the charm content of the photon. The remaining ln(p^/m 2 ) terms connected with the 
collinear singularities at the proton leg are higher order in a s and are definitely much less important. 

In Fig.lc we have plotted also the sum of the direct and the resolved components for the massless 
theory including NLO corrections (denoted NLO massless). The higher order corrections to the 
resolved process increase the theoretical prediction by 35% compared to the NLO/LO massless 
result. Following the discussion above this contribution should correspond to a large extent to the 
NNLO direct cross section in the massive scheme which has not been calculated yet. 

In Fig. 2 the differential cross sections d 2 <r/dydp 2 L are shown for p± = 10 GeV as a function of 
the rapidity y for the direct (Fig. 2a) and the resolved contributions (Fig. 2b). Comparing the results 
in the massive and the massless approaches we observe the same pattern as for the p± distribution, 
Figs. la and b. From Fig. 2a we can infer that both the LO and NLO predictions for the direct 
contribution in the massive and massless approaches are very close together. In Fig. 2b, where the 
resolved y differential cross sections are plotted we see the very much reduced cross section in the 
massive theory and the appreciable cross section in the massless approach, where the NLO prediction 
is in general larger than the LO result. 

Fig. 2c finally shows the sum of direct and resolved contributions in the massless and massive 
schemes. The same relation between massive and massless schemes is seen here as in the p± distri- 
bution in Fig.lc. 
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3. Conclusion 



We have observed that the cross section in the massless approach is approximately 110% larger 
than in the massive scheme. The main part of that deviation can be attributed to the difference 
found already in the comparison of the NLO/LO massless cross section with the NLO/LO massive 
result. This difference originates from several sources which have been disentangled. The remaining 
increase of the NLO massless prediction is due to the NLO corrections to the resolved process. Such 
a contribution corresponds mainly to a NNLO correction to the direct part which has not been 
calculated yet. In this sense the massless theory includes a significant part of the NNLO correction 
and should therefore yield a more reliable prediction than the massive scheme. This conclusion is 
supported by the observation that the massless prediction is much more stable under scale variation 
than the massive result. 



Compared to a similar study for b production in pp processes fig] we find much larger deviations 
between the massive and the massless approaches. This can be attributed to the much stronger 
influence of charm in the photon as compared to the proton caused by the point-like component of 
the photon structure function. Consequently a measurement of large-pj_ charm production at HERA 
might provide information about the charm content of the photon (cf. |21| ). Such measurements 
will be very instructive since theoretical opinions on that issue are rather divided p4| ]. 

Our study was based on a scale-independent 5-type fragmentation function of the charm quarks. 
In the future we shall extend this by using more realistic forms for the fragmentation including 
evolution to higher scales. This latter effect in the fragmentation might compensate part of the 
ln(p^/m 2 ) absorption in the massless theory and might thus bring the massless prediction nearer 
to the massive result. 

Note added: After completion of our work we have been informed about a similar analysis by Cacciari 
and Greco [25|. 
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Fig. la: Transverse momentum distributions of ep — > c/c + X averaged over the rapidity range 
\y\ < 1.5 for the direct contribution in the massless and massive schemes. 
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Fig. lb: Transverse momentum distributions of ep — ► c/c + X averaged over the rapidity range 
\y\ < 1.5 for the resolved contribution in the massless and massive schemes. 
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Fig.lc: Transverse momentum distributions of ep — > c/c + X averaged over the rapidity range 
\y\ < 1.5 in the massless and massive schemes. Shown is the sum of the NLO massive direct and the 
LO massive resolved contributions (NLO/LO massive), the sum of the NLO massless direct and the 
LO massless resolved contributions (NLO/LO massless), and the sum of the NLO massless direct 
and the NLO massless resolved contributions (NLO massless). 
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Fig.2a: Rapidity distributions of ep — > c/c + X at p± = 10 GeV for the direct contribution in the 
massless and massive schemes. 
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Fig.2b: Rapidity distributions of ep — > c/c + X at p± = 10 GeV for the resolved contribution in 
the massless and massive schemes. 
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Fig.2c: Rapidity distributions of ep — > c/c + X af = 10 GeV for the resolved contribution in the 
massless and massive schemes. Shown is the sum of the NLO massive direct and the LO massive 
resolved contributions (NLO/LO massive), the sum of the NLO massless direct and the LO massless 
resolved contributions (NLO/LO massless), and the sum of the NLO massless direct and the NLO 
massless resolved contributions (NLO massless). 
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